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Cosmological Structure Formation 


1.1 Introduction 


Cosmology has finally become a mature science during the past decade, with 
predictions now routinely confirmed by observations. The modern cosmolog¬ 
ical theory is known as ACDM - CDM for Cold Dark Matter, particles that 
moved sluggishly in the early universe and thereby preserved fluctuations 


down to small scales (|Blumenthal et al. 

1984 

see Fig. 1), and A for the cos- 

mological constant (e.g., 

Lahav et al. 

1991 

. A wide variety of large-scale 


astronomical observations - including the Cosmic Microwave Background 
(CMB), measurements of baryon acoustic oscillations (BAO), gravitational 
leasing, the large-scale distribution of galaxies, and the properties of galaxy 
clusters - agree very well with the predictions of the ACDM cosmology. 


Like the standard model of particle physics, the ACDM standard cosmo¬ 
logical model requires the determination of a number of relevant cosmolog¬ 
ical parameters, and the theory does not attempt to explain why they have 
the measured values - or to explain the fundamental nature of the dark mat¬ 
ter and dark energy. These remain challenges for the future. But the good 
news is that the key cosmological parameters are now all determined with 
unprecedented accuracy, and the six-parameter ACDM theory provides a 
very good match to all the observational data including the 2015 Planck 


temperature and polarization data (Planck Collaboration et al. 2015a) 


Within uncertainties less than 1%, the Universe has critical cosmic den¬ 
sity - i.e., Dtotai = 1-00 and the Universe is Euclidean (or “flat”) on large 
scales. The expansion rate of the Universe is measured by the Hubble pa¬ 
rameter h = 0.6774 ± 0.0046, and ^matter = 0.3089 ± 0.0062; this leads to 
the age of the Universe to = 13.799 ± 0.021 Gyr. The power spectrum nor¬ 
malization parameter is cg = 0.816 ± 0.009, and the primordial fluctuations 
are consistent with a purely adiabatic spectrum of fluctuations with a spec- 
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tral tilt Ug = 0.968 ± 0.006, as predicted by single-field inflationary models 
(Planck Collaboration et al., 2015a). The same cosmological parameters that 
are such a good match to the CMB observations also predict the observed 
distribution of density fluctuations from small scales probed by the Lyman 
alpha fores10to the entire horizon, as shown in Fig. 2. The near-power-law 
galaxy-galaxy correlation function at low redshifts is now known to be a 
cosmic coincidence (Watson et al., j2011 ). I was personally particularly im¬ 
pressed that the evolution of the galaxy-galaxy correlations with redshift 
predicted by ACDM (Kravtsov et ah, 2004) turned out to be in excellent 
agreement with the subsequent observations (e.g., Conroy et ah, 2006). 

For non-astronomers, there should be a more friendly name than ACDM 
for the standard modern cosmology. Since about 95% of the cosmic density 
is dark energy (either a cosmological constant with Da = 0.69 or some 
dynamical field that plays a similar cosmic role) and cold dark matter with 
hlcDM = 0.26, I recommend the simple name “Double Dark Theory” for the 
modern cosmological standard model (jPrimack and Abrams 2006 Abrams 


and Primack, 2011). The contribution of ordinary baryonic matter is only 


Db = 0.05. Only about 10% of the baryonic matter is in the form of stars 
or gas clouds that emit electromagnetic radiation, and the contribution of 
what astronomers call “metals” - chemical elements heavier than helium ~ 
to the cosmic density is only Dmetais ~ 0.0005, most of which is in white 
dwarfs and neutron stars (Fukugita and Peebles, 2004). The contribution 
of neutrino mass to the cosmic density is 0.002 < Dj/ < 0.005, far greater 
than Djnetais • Thus our bodies and our planet are made of the rarest form of 
matter in the universe: elements forged in stars and stellar explosions. 

Potential challenges to ACDM on large scales come from the tails of the 
predicted distribution functions, such as CMB cold spots and massive clus¬ 
ters at high redshifts. However, the existing observations appear to be consis¬ 
tent thus far with predictions of standard ACDM with standard primordial 


power spectra; non-Gaussian initial conditions are not required (Planck Col¬ 


laboration et ah, 2015b). Larger surveys now underway may provide more 


stringent tests. 


1.2 Large-Scale Structure 

Large, high-resolution simulations permit detailed predictions of the distri¬ 
bution and properties of galaxies and clusters. From 2005-2010, the bench¬ 


mark simulations were Millennium-I (Springel et ah, 2005) and Millennium' 


^ The Lya forest is the many absorption lines in quasar spectra due to clouds of neutral 
hydrogen along the line of sight to the quasar. 
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Figure 1.1 The origin of the CDM spectrum of density fluctuations. Left 
panel: Fluctuations corresponding to mass scales 10® M©, 10® Mq, etc., 
grow proportionally to the square of scale factor a when they are outside 
the horizon, and when they enter the horizon (cross the horizontal dashed 
line) the growth of the fluctuation amplitude S is much slower if they enter 
when the Universe is radiation dominated (i.e., a < Oeq). Fluctuations on 
mass scales > 10^^Mq enter the horizon after it becomes matter dominated, 
so their growth is proportional to scale factor a; that explains the larger 
separation between amplitud es for such higher-mass fluctuat ions. (From a 
1983 con f erence presentation Primack and Blumenthal| ( 1984 ), reprinted in 
Primackj (1984|.) Right panel: The resulting CUM fluctuation spectrum 
= AM/M) is contrasted with a 5 oc k® white noise (Poisson) 
spectrum with the same power at all wavelengths (where wave number 
K is as usual related to wavelength A by k = 27r/A), and with the hot 
dark matter spectrum if the dark matter were light neutrinos (ly) which 


i s cut off on galaxy mass scales by free-steaming. (From Blumenthal et al. 
(1984). This calcnlation assumed that the primordial fluctuations are scaie- 
invariant (Zel’dovich) and that Umatter = 1 and Hnbble parameter h = 1.) 


II (Boylan-Kolchin et al., 2009), which have been the basis for more than 


400 papers. However, these simulations used first-year Wilkinson Microwave 
Anisotropy Probe (WMAP) cosmological parameters, including as = 0.90, 
that are now in serious disagreement with the latest observations. Improved 
cosmological parameters, simulation codes, and computer power have per¬ 
mitted more accurate simulations (]Kuhlen et al. 2012[ Skillman et al.| 2014) 


including Bolshoi (Klypin et al. 


et al. 

2012 

Riebe et al. 

2013 


Planck cosmological parameters (Klypin et al. 


2011) and BigBolshoi/MultiDark (Prada 
cently be 

2014b[ )P1 


Dark matter halos can be characterized in a number of ways. A common 
one is by mass, but the mass attributed to a halo depends on a number of 


^ The web address for the MultiDark simulation data center is 
http://www.cosmosim.org/cms/simulations/multidark-project/ 
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Mass scale M [Msolar] 


Figure 1.2 The r.m.s. mass variance AM/M predicted by ACDM com¬ 
pared with observations, from CMB and the Atacama Cosmology Telescope 
(ACT) on large scales, brightest cluster galaxy weak lensing, clusters, the 
SDSS galaxy distribution, to the Lyman alpha forest on small scales. This 
figure highlights the consistency of power spectrum measurements by an 
array of cosmological probes over a large range of scales. (Redrawn from 


Fig. 5 in Hlozek et al. (2012), which gives the sources of the data.) 


factors including how the outer edge of the halo is defined; popular choices 
include the spherical radius within which the average density is either 200 
times critical density or the virial density (which depend on redshift). Prop¬ 
erties of all the halos in many stored time steps of both the Bolshoi and 
BigBolshoi/MultiDark simulations are available on the web in the Multi- 
Dark database.^ For many purposes it is more useful to characterize halos 
by their maximum circular velocity Ihnax; which is defined as the maximum 
value of [GM{< r)/r]^/^, where G is Newton’s constant and M(< r) is 
the mass enclosed within radius r. The reason this is useful is that Imax is 
reached at a relatively low radius rmaxj closer to the central region of a halo 
where stars or gas can be used to trace the velocity of the halo, while most 
of the halo mass is at larger radii. Moreover, the measured internal veloc- 
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ity of a galaxy (line of sight velocity dispersion for early-type galaxies and 
rotation velocity for late-type galaxies) is closely related to its luminosity 
according to the Faber-Jackson and Tully-Fisher relations. In addition, after 
a subhalo has been accreted by a larger halo, tidal stripping of its outer parts 
can drastically reduce the halo mass but typically decreases I4iax much less. 
(Since the stellar content of a subhalo is thought to be determined before it 
was accreted, some authors define Imax to be the peak value at any redshift 
for the main progenitor of a halo.) Because of the observational connec¬ 
tion between larger halo internal velocity and brighter galaxy luminosity, a 
common simple method of assigning galaxies to dark matter halos and sub¬ 
halos is to rank order the galaxies by luminosity and the halos by Vmaxj and 


then match them such that the number densities are comparable (Kravtsov 


et al., 2004, Tasitsiomi et al.[ 2004; Conroy et al., 2006). This is called “halo 
abundance matching” or (more modestly) “sub-halo abundance matching” 


(SHAM) (Reddick et al., 2014). Halo abundance matching using the Bolshoi 


simulation predicts galaxy-galaxy correlations (which are essentially counts 
of the numbers of pairs of galaxies at different separation distances) that are 
in good agreement with the Sloan Digital Sky Survey (SDSS) observations 
( ]Trujillo-Gomez et al. 2011; Reddick et al. 


2013). 


Abundance matching with the Bolshoi simulation also predicts galaxy 


velocity-mass scaling relations consistent with observations (Trujillo-Gomez 


et al. 2011), and a galaxy velocity function in good agreement with ob¬ 


servations for maximum circular velocities l^aax ~ 100 km/s, but higher 
than the HI Parkes All Sky Survey (HIPASS) and the Arecibo Legacy Fast 


ALFA (ALFALFA) Survey radio observations (Zwaan et al., 2010 Papaster- 


gis et al.] 2011) by about a factor of 2 at 80 km/s and a factor of 10 at 50 


km/s. This either means that these radio surveys are increasingly incom¬ 
plete at lower velocities, or else AGDM is in trouble because it predicts far 
more small-14iax halos than there are observed low-P galaxies. A deeper 
optical survey out to 10 Mpc found no disagreement between 14iax predic¬ 
tions and observations for 14 iax > 60 km/s, and only a factor of 2 excess 
of halos compared to galaxies at 40 km/s (Klypin et al., 2014a). This may 
indicate that there is no serious inconsistency with theory, since for K 30 
km/s reionization and feedback can plausibly explain why there are fewer 


observed galaxies than dark matter halos (Bullock et al., 2000; Somerville 


20021 Benson et al. 

2003 

Kravtsov, 2010 Wadepuhl and Springell 2011 

Sawala et al.l 2012|, and also the observed scaling of metallicity with galaxy 

mass ( 

Dekel and Wo( 

3, 2003 

Woo et al., 2008; Kirby et al. 

2011). 


The radial dark matter density distribution in halos can be approximately 
fit by the simple formula pnfw = ^PsX~^{1 -\-x)~‘^, where x = r/r^ (Navarro 
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et al., 1996), and the “concentration” of a dark matter halo is defined as 
.g where i?vir is the virial radius of the halo. When we first under- 


C — RtjIy/ Rg 

stood that dark matter halos form with relatively low concentration C ~ 4 
and evolve to higher concentration, we suggested that “red” galaxies that 
shine mostly by the light of red giant stars because they have stopped form¬ 
ing stars should be found in high-concentration halos while “blue” galaxies 
that are still forming stars should be found in younger low-concentration 


halos (Bullock et ah, 2001). This idea was recently rediscovered by Hearin 


and Watson, who used the Bolshoi simulation to show that this leads to re¬ 
markably accurate predictions for the correlation functions of red and blue 


galaxies (Hearin and Watson, 2013 Hearin et al. 2014). 


The Milky Way has two rather bright satellite galaxies, the Large and 
Small Magellanic Clouds. It is possible using sub-halo abundance matching 
with the Bolshoi simulation to determine the number of Milky-Way-mass 
dark matter halos that have subhalos with high enough circular velocity 
to host such satellites. It turns out that about 55% have no such subha¬ 
los, about 28% have one, about 11% have two, and so on (]Busha et al 


2011a). Remarkably, these predictions are in excellent agreement with an 


analysis of observations by the Sloan Digital Sky Survey (SDSS) (Liu et al. 


2011). The distribution of the relative velocities of central and bright satel¬ 


lite galaxies from SDSS spectroscopic observations is also in very good 


agreement with the predictions of the Millennium-H simulation (Tollerud 


et al.| 2011), and the Milky Way’s lower-luminosity satellite population 


is not unusual (Strigari and Wechsler 

2012). Considered in a cosmologi- 

cal context, the Magellanic cloud 
about the last Gyr (Besla et al. 

s are 

[m2 

likely to have been accreted within 
, and the Milky Way halo mass is 

1.2lQ;4(stat.)±0.3(sys.)xl0^^MQ 

Busha et ah, 2011b). 


1.3 Galaxy Formation 

At early times, for example the CMB epoch about 400,000 years after the big 
bang, or on very large scales at later times, linear calculations starting from 
the ACDM fluctuation spectrum allow accurate predictions. But on scales 
where structure forms, the fluctuations have grown large enough that they 
are strongly nonlinear, and we must resort to simulations. The basic idea is 
that regions that start out with slightly higher than average density expand 
a little more slowly than average because of gravity, and regions that start 
out with slightly lower density expand a little faster. Nonlinear structure 
forms by the process known by the somewhat misleading name “gravitational 
collapse” - misleading because what really happens is that when positive 
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rrlilkyWay 
100,000 Light Years 



Milky Way Dark Matter Halo 
1,500,000 Light Years 



Figure 1.3 The stellar disk of a large spiral galaxy like the Milky Way is 
about 100,000 light years across, which is tiny compared with the dark 
matter halo of such a galaxy (from the Aquarius dark matter simulation 
Springel et al. 2008|), and even much smaller compared with the large-scale 
cosmic web (from the Bolshoi simulation Klypin et al. 20111. 


fluctuations have grown sufficiently that they are about twice as dense as 
typical regions their size, they stop expanding while the surrounding universe 
keeps expanding around them. The result is that regions that collapse earlier 
are denser than those that collapse later; thus galaxy dark matter halos are 
denser than cluster halos. The visible galaxies form because the ordinary 
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baryonic matter can radiate away its kinetic energy and fall toward the 
centers of the dark matter halos; when the ordinary matter becomes dense 
enough it forms stars. Thus visible galaxies are much smaller than their 
host dark matter halos, which in turn are much smaller than the large scale 
structure of the cosmic web, as shown in Fig. 3. 

Astronomical observations represent snapshots of moments long ago when 
the light we now observe left distant astronomical objects. It is the role 
of astrophysical theory to produce movies - both metaphorical and actual 
- that link these snapshots together into a coherent physical picture. To 
predict cosmological large-scale structures, it has been sufficient to treat all 
the mass as dark matter in order to calculate the growth of structure and 
dark matter halo properties. But hydrodynamic simulations - i.e., including 
baryonic matter - are necessary to treat the formation and evolution of 
galaxies. 
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Figure 1.4 The Star-Forming Band on a diagram of baryon density nt ver¬ 
sus the three-dimensional r.m.s. velocity dispersion V and virial tempera¬ 
ture T for structures of various sizes in the universe, where T = /rF^/3/c, 
p. is mean molecular weight (« 0.6 for ionized primordial H -|- He) and 
k is Boltzmann’s constant. Below the No Metals and Solar Metals cooling 
curves, the cooling timescale is more rapid than the gravitational timescale. 
Dots are groups and clusters. Diagonal lines show the halo masses in units 
of Mq. (This is Fig. 3 in Blumenthal et ah, 1984, with the Star-Forming 
Band added.) 
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An old criticism of ACDM has been that the order of cosmogony is wrong: 
halos grow from small to large by accretion in a hierarchical formation the¬ 
ory like ACDM, but the oldest stellar populations are found in the most 
massive galaxies - suggesting that these massive galaxies form earliest, a 


phenomenon known as “downsizing” (Cowie et ah, 1996). The key to ex¬ 


plaining the downsizing phenomenon is the realization that star formation 
is most efficient in dark matter halos with masses in the band between about 


10^^ and IO^^Mq (Eig. 1 bottom in 

Behroozi et al. 

2011 

>). This goes back 

at least as far as the original Cold Dark Matter paper ( 

llumenthal et al. 


1984): see Fig. 4. A dark matter halo that has the total mass of a cluster 


of galaxies today will have crossed this star-forming mass band at an early 
epoch, and it will therefore contain galaxies whose stars formed early. These 
galaxies will be red and dead today. A less massive dark matter halo that 
is now entering the star-forming band today will just be forming significant 
numbers of stars, and it will be blue today. The details of the origin of the 
star-forming band are still being worked out. Back in 1984, we argued that 
cooling would be inefficient for masses greater than about IO^^Mq because 
the density would be too low, and inefficient for masses less than about 
10®Mq because the gas would not be heated enough by falling into these 


small potential wells. Now we know that reionization, supernovae (Dekel and 


Silk 1986), and other energy input additionally impedes star formation for 
halo masses below about IO^^Mq, and feedback from active galactic nuclei 
(AGN) additionally impedes star formation for halo masses above about 

Early simulations of disk galaxy formation found that the stellar disks 


had much lower rotation velocities than observed galaxies (Navarro and 


Steinmetz, 2000). This problem seemed so serious that it became known 
as the “angular momentum catastrophe.” A major cause of this was ex¬ 
cessive cooling of the gas in small halos before they merged to form larger 
galaxies (Mailer and Dekel, 2002). Simulations with better resolution and 
more physical treatment of feedback from star formation appear to resolve 


this problem. In particular, the Eris cosmological simulation (Guedes et al. 


2011) produced a very realistic spiral galaxy, as have many simulations since 


then. Somerville and Dave (2014) is an excellent recent review of progress in 
understanding galaxy formation. In the following I summarize some of the 
latest developments. There are now two leading approaches to simulating 
galaxies: 


• Low resolution, ~ 1 kiloparsecs. The advantages of this approach are 
that it is possible to simulate many galaxies and study galaxy popula- 
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tions and their interactions with the circumgalactic and intergalactic me¬ 
dia, but the disadvantages are that we learn relatively little about how 
galaxies themselves form and evolve at high redshifts. The prime exam¬ 


ples of this approach now are the Illustris (Vogelsberger et ah, 2014b) 


and EAGLE ( Schaye et al.[ 2015[ ) simulations. Like semi-analytic models 
of galaxy formation (reviewed in Benson, 2010), these projects adjusted 
the parameters governing star-formation and feedback processes in or¬ 
der to reproduce key properties of galaxies at the present epoch, redshift 
z = 0. The Illustris simulation in a volume (106.5Mpc)^ forms ~ 40,000 
galaxies at the present epoch with a reasonable mix of elliptical and spi¬ 


ral galaxies that have realistic appearances (Snyder et al. 2015), obey 


observed scaling relations, have the observed numbers of galaxies as a 
function of their luminosity, and were formed with the observed cosmic 


star formation rate (Vogelsberger et al., 2014a). It forms massive compact 


galaxies by redshift z = 2 via central starbursts in major mergers of gas- 


rich galaxies or else by assembly at very early times (Wellons et al., 2015). 


Remarkably, the Illustris simulation also predicts a population of damped 
Lyman a. absorbers (DLAs, small-galaxy-size clouds of neutral hydrogen) 
that agrees with some of the key observational properties of DLAs ([Bird 


et ah, 2014, 2015). The EAGLE simulation in volumes up to (lOOMpc)^ 


reproduces the observed galaxy mass function from 10® to IO^^Mq at a 
level of agreement close to that attained by semi-analytic models (|Schaye 


et ah 


2015), and the observed atomic and molecular hydrogen content of 
3 (iRahmati et ah I 


2015; Lagos et ah, 2015). 


galaxies out to z 

High resolution, ~ 10s of parsecs. The advantages are that it is possible 
to compare simulation outputs in detail with high-redshift galaxy images 
and spectra to discover the drivers of morphological changes as galaxies 
evolve, but the disadvantage is that such simulations are so expensive 
computationally that it is as yet impossible to run enough cases to create 
statistical samples. Leading examples of this approach are FIRE simula¬ 


tions led by Phil Hopkins (e.g., Hopkins et ah, 2014) and the ART simula¬ 


tion suite led by Avishai Dekel and myself (e.g., Zolotov et ah, 2014). We 


try to compensate for the small number of high-resolution simulations by 
using simulation outputs to tune semi-analytic models, which in turn use 
cosmological dark-matter-only simulations like Bolshoi to follow the evo¬ 
lution of ~ 10® galaxies in their cosmological context (e.g. 

2014b|a Brennan et al 


Porter et ah 


2015). 


The high-resolution FIRE simulations, based on the GIZMO smooth par¬ 
ticle hydrodynamics code ( Hopkinsf 2014) with supernova and stellar feed- 
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back, including radiative feedback (RF) pressure from massive stars, treated 
with zero adjusted parameters, reproduce the observed relation between stel¬ 
lar and halo mass up to Mhaio ~ IO^^Mq and the observed star formation 
rates Hopkins et al. (2014). FIRE simulations predict covering fractions of 
neutral hydrogen with column densities from 10^^cm“^ (Lyman limit sys¬ 
tems, LLS) to > 10^*^'^cm“^ (DLAs) in agreement with observations at red- 
shifts z=2-2.5 (Faucher-Giguere et ah, 2015); this success is a consequence 
of the simulated galactic winds. FIRE simulations also correctly predict the 


observed evolution of the decrease of metallicity with stellar mass (Ma et al. 


2015), and produce dwarf galaxies that appear to agree with observations 
( jOhorbe et al. 2015) as we will discuss in more detail below. 

The high-resolution simulation suite based on the ART adaptive mesh 


refinement (AMR) approach (Kravtsov et ah, 1997; Ceverino and Klypin 


2009) incorporates at the sub-grid level many of the physical processes rel¬ 


evant for galaxy formation. Our initial group of 30 zoom-in simulations of 
galaxies in dark matter halos of mass (1 — 30) x IO^^Mq at redshift z = 1 
were run at 35-70 pc maximum (physical) resolution ([Ceverino et al 


2012 


2015a). The second group of 35 simulations (VELAOl to VELA35) with 17.5 
to 35 pc resolution of halos of mass (2 — 20) x IO^^Mq at redshift z = lhave 
now been run three times with varying inclusion of radiative pressure feed- 


back (none, UV, UV-I-IR), as described in 

Ceverino et ah ( 

2014 

). RF pressure 

including the effects of stellar winds ( 

Hopkins et ah 

2012 

2014) captures es- 


sential features of star formation in our simulations. In particular, RF begins 
to affect the star-forming region as soon as massive stars form, long before 
the first supernovae occur, and the amount of energy released in RF greatly 


exceeds that released by supernovae (Ceverino et ah, 2014; Trujillo-Gomez 


et ah , 2015). In addition to radiation pressure, the local UV flux from young 


star clusters also affects the cooling and heating processes in star-forming 
regions through photoheating and photoionization. We use our Sunrise code 
dJonsson 2006 Jonsson et ah, 2006, 2010 Jonsson and Primack, 2010) to 


make realistic images and spectra of these simulated galaxies in many wave¬ 
bands and at many times during their evolution, including the effects of stel¬ 
lar evolution and of dust scattering, absorption, and re-emission, to compare 
with the imaging and photometry from CANDEL^and other surveys - see 
Fig. 5 for examples including the effect of CANDELization (reducing the 
resolution and adding noise) to allow direct comparison with HST images. 

In comparing our simulations with HST observations, especially those 
from the CANDELS and 3D-HST surveys, we are hnding that the simula- 


3 


CANDELS, the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, was the 
largest-ever Hubble Space Telescope survey, see http://caiidels.ucolick.org/ 
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V26 zs3.6 dSF 

V26 zs2.7 cSF, ex situ clump 

V26 zb2.3 cSF, in situ clumps 

r 

V26zb1.2 cQ, minor merger 






Figure 1.5 Face-on images of Vela26 simulated galaxy with UV radiation 
pressure feedback, at four redshifts (a) z = 3.6 when it is diffuse and star 
forming (dSF); (b) z = 2.7 when it has become compact and star forming 
(cSF) with a red ex-situ clump; (c) z = 2.3 still cSF, now with in situ 
clumps apparent in the V-band image; (d) compact and quenched (cQ) 
during a minor merger, with tidal features visible in the V-band image. 
Top panels: three-color composite images at high resolution; bottom pan¬ 
els: CANDELized V and H band images. The observed V band images 
correspond to ultraviolet radiation from massive young stars in the galaxy 
rest frame, while the observed H band images show optical light from the 
entire stellar population including old stars. The CANDELS survey took 
advantage of the infrared capability of the Wide Field Camera 3, installed 
on the last service visit to HST in 2009. 


tions can help us interpret a variety of observed phenomena that we now 
realize are important in galaxy evolution. One is the formation of com¬ 


pact galaxies. Analysis of CANDELS images suggested (Barro et al., 2013 


2014a|b) that diffuse star-forming galaxies become compact galaxies (blue 


nuggets) which subsequently quench (red nuggets). We see very similar be¬ 


havior in our VELA simulations with UV radiative feedback (Zolotov et al. 


2014, see Eigure 2), and we have identified in our simulations several mecha¬ 


nisms that lead to compaction often followed by rapid quenching, including 
major gas-rich mergers, disk instabilities often triggered by minor mergers, 


and opposing gas flows into the central galaxy (Danovich et ah, 2015). 


Another aspect of galaxy formation seen in HST observations is massive 

and references 


star-forming clumps (Guo et ah, 2012 Wuyts et ah 2013 


therein), which occur in a large fraction of star-forming galaxies at redshifts 


z = 1-3 (Guo et ah, 2015). In our simulations there are two types of clumps. 


Some are a stage of minor mergers - we call those ex situ clumps. A majority 
of the clumps originate in situ from violent disk instabilities (VDI) in gas- 


rich galaxies (Ceverino et ah, 2012 Moody et ah, 2014 Mandelker et al 


2014). Some of these in situ clumps are associated with gas instabilities that 


help to create compact spheroids, and some form after the central spheroid 
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and are associated with the formation of surrounding disks. We find that 
there is not a clear separation between these processes, since minor mergers 


2014). 


often trigger disk instabilities in our simulations (Zolotov et al. 

Star-forming galaxies with stellar masses M* 3 x 10 ®Mq at z > 1 have 
recently been shown to have mostly elongated (prolate) stellar distri¬ 


butions (van der Wei et al., 2014) rather than disks or spheroids, based on 


their observed axis ratio distribution. In our simulations this occurs because 


most dark matter halos are prolate especially at small radii (Allgood et al. 


2006), and the first stars form in these elongated inner halos; at lower red- 


shifts, as the stars begin to dominate the dark matter, the galaxy centers 


become disky or spheroidal (Ceverino et al., 2015b). 


Both the FIRE and ART simulation groups and many others are partici¬ 
pating in the Assembling Galaxies of Resolved Anatomy (AGORA) collabo¬ 
ration (Kim et al., 2014) to run high-resolution simulations of the same initial 
conditions with halos of masses 10^®, 10^^, 10^^, and IO^^Mq at z = 0 with 
as much as possible the same astrophysical assumptions. AGORA cosmo¬ 
logical runs using different simulation codes will be systematically compared 
with each other using a common analysis toolkit and validated against ob¬ 
servations to verify that the solutions are robust - i.e., that the astrophysical 
assumptions are responsible for any success, rather than artifacts of partic¬ 
ular implementations. The goals of the AGORA project are, broadly speak¬ 
ing, to raise the realism and predictive power of galaxy simulations and the 
understanding of the feedback processes that regulate galaxy “metabolism.” 

It still remains to be seen whether the entire population of galaxies can 
be explained in the context of ACDM. A concern regarding disk galaxies is 
whether the formation of bulges by both galaxy mergers and secular evo¬ 
lution will prevent the formation of as many pure disk galaxies as we see 
in the nearby universe (Kormendy and Fisher, 2008). A concern regarding 
massive galaxies is whether theory can naturally account for the relatively 
large number of ultra-luminous infrared galaxies. The bright sub-millimeter 
galaxies were the greatest discrepancy between our semi-analytic model pre¬ 


dictions compared with observations out to high redshift (Somerville et al. 


2012). This could possibly be explained by a top-heavy stellar initial mass 


function, or perhaps more plausibly by more realistic simulations including 


self-consistent treatment of dust (Hayward et al., 2011, 2013). Glearly, there 


is much still to be done, both observationally and theoretically. It is possi¬ 
ble that all the potential discrepancies between AGDM and observations of 
relatively massive galaxies will be resolved by better understanding of the 
complex astrophysics of their formation and evolution. But small galaxies 
might provide more stringent tests of AGDM. 
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1.4 Smaller Scale Issues: Cusps 


Cusps were perhaps the first potential discrepancy pointed out between the 
dark matter halos predicted by CDM and the observations of small galaxies 


that appeared to be dominated by dark matter nearly to their centers (Flores 


and Primack 

1994 

Moore 

1994 




that the central density of dark matter halos behaves roughly as p ~ r 
As mentioned above, dark matter halos have a density distribution that 
can be roughly approximated as pnfw = 4p<jX“^(l + where x = r/rs 
(Navarro et al., 1996). But this predicted r~^ central cusp in the dark matter 


distribution seemed inconsistent with published observations of the rotation 
velocity of neutral hydrogen as a function of radius. 


In small galaxies with significant stellar populations, simulations show 
that central starbursts can naturally produce relatively fiat density profiles 


(Covernato et al. 

2010 

201 ^ i 

^ontzen and Governato 

2012 

Teyssier et al. 

2013t Brooks, 2014t Brooks and Zolotov 2014| Madau et al. 

2014| Ohorbe 

et al., 2015, Nipoti and Binney 

2015). Gas cools into the galaxy center and 


becomes gravitationally dominant, adiabatically pulling in some of the dark 


matter (Blumenthal et ah, 1986; Gnedin et al., 2011). But then the gas is 


driven out very rapidly by supernovae and the entire central region expands, 
with the density correspondingly dropping. Several such episodes can occur, 
producing a more or less constant central density consistent with observa¬ 
tions, as shown in Fig. 6. The figure shows that galaxies in the THINGS 
sample are consistent with ACDM hydrodynamic simulations. But simu¬ 
lated galaxies with stellar mass less than about 3 x 10 ®Mq may have cusps, 
although Ohorbe et al. (2015) found that stellar effects can soften the cusp 
in even lower-mass galaxies if the star formation is extended in time. The 
observational situation is unclear. In Sculptor and Fornax, the brightest 
dwarf spheroidal satellite galaxies of the Milky Way, stellar motions may 
imply a flatter central dark matter radial profile than p ~ r“^(Walker and 


Peharrubia 

2011 

Amorisco and Evans 

2012 

|Jardel and Gebhardt[ 

2012D. 

However, other papers have questioned this (jJardel and Gebhardt 

2013 


Breddels and Helmi, 2013, 2014; Richardson and Fairbairn 


2014). 


Will baryonic effects explain the radial density distributions in larger low 
surface brightness (LSB) galaxies? These are among the most common galax¬ 
ies. They have a range of masses but many have fairly large rotation veloc¬ 
ities indicating fairly deep potential wells, and many of them may not have 
enough stars for the scenario just described to explain the observed rota- 
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Figure 1.6 Dark matter cores are generated by baryonic effects in galaxies 
with sufficient stellar mass. The slope ex of the dark matter central den¬ 
sity profile r“ is plotted vs. stellar mass measured at 500 parsecs from 


simulations described in 

Pontzen and Governato 

(I 2 OI 2 I). The solid NFW 

curve assumes the halo concentrations given by (JV 

laccio et al. 

2007 

). Large 


crosses: halos with > 5 x 10® dark matter particles; small crosses: > 5 x 10^ 
particles. Squares represent galaxies observed by The HI Nearby Galaxy 
Survey (THINGS). (Fig. 3 in Pontzen and Governato 2014.) 


2002) and the observed diversity of rotation curves (Maccio et ah, 2012b 


Oman et al. 2015)? This is a challenge for galaxy simulators. 


Some authors have proposed that warm dark matter (WDM), with initial 
velocities large enough to prevent formation of small dark matter halos, 
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could solve some of these problems. However, that does not appear to work: 
the systematics of galactic radial density profiles predicted by WDM do not 


at all match the observed ones (Kuzio de Naray et ah, 2010, Maccio et al 


2012a, 2013). WDM that’s warm enough to affect galaxy centers may not 


permit early enough galaxy formation to reionize the universe (Governato 


et al., 2015). Yet another constraint on WDM is the evidence for a great deal 


of dark matter substructure in galaxy halos (Zentner and Bullock, 2003), 
discussed further below. 


1.5 Smaller Scale Issues: Satellite Galaxies 


As the top panel of Fig. 3 shows, ACDM predicts that there are many fairly 
massive subhalos within dark matter halos of galaxies like the Milky Way 
and the Andromeda galaxy, more than there are observed satellite galaxies 


Klypin et ah 

1999 

Moore et ah 

1999 


1999). This is not obviously a problem 


for the theory since reionization, stellar feedback, and other phenomena are 
likely to suppress gas content and star formation in low-mass satellites. As 
more faint satellite galaxies have been discovered, especially using multicolor 
information from SDSS observations, the discrepancy between the predicted 
and observed satellite population has been alleviated. Many additional satel¬ 
lite galaxies are predicted to be discovered by deeper surveys (e.g., Bullock 


et ah, 2010), including those in the Southern Hemisphere seen by the Dark 


Energy Survey (The DES Collaboration et ah, 2015) and eventually the 


Large Synoptic Survey Telescope (LSST). 

However, a potential discrepancy between theory and observations is the 
“too big to fail” (TBTF) problem (Boylan-Kolchin et ah, 2011, 2012). The 


Via Lactea-H high-resolution dark-matter-only simulation of a Milky Way 


size halo ( 

Diemand et ah 

2007 

2008) and the six similar Aquarius Simula- 

tions ( 

Springel et ah 

2008) all have several subhalos that are too dense in 


their centers to host any observed Milky Way satellite galaxy. The brightest 
observed dwarf spheroidal (dSph)satellites all have 12 Fmax ~ 25 km/s. 
But the Aquarius simulations predict at least 10 subhalos with Vnax > 25 
km/s. These halos are also among the most massive at early times, and thus 
are not expected to have had their star formation greatly suppressed by 
reionization. They thus appear to be too big to fail to become observable 


satellites (Boylan-Kolchin et ah, 2012). 


The TBTF problem is closely related to the cusp-core issue, since TBTF is 
alleviated by any process that lowers the central density and thus the internal 
velocity of satellite galaxies. Many of the papers finding that baryonic effects 
remove central cusps cited in the previous section are thus also arguments 
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against TBTF. A recent simulation of regions like the Local Group found 
the number, internal velocities, and distribution of the satellite galaxies to 


be very comparable with observations (Sawala et ah, 2014). 


Perhaps there is additional physics beyond ACDM that comes into play on 
small scales. One possibility that has been investigated is warm dark matter 
(WDM). A simulation like Aquarius but with WDM has fewer high-lmax 
halos (Lovell et ah, 2012). But it is not clear that such WDM simulations 
with the lowest WDM particle mass allowed by the Lyman alpha forest 


and other observations (Viel et ah, 2013 Horiuchi et ah 2014) will have 


enough substructure to account for the observed faint satellite galaxies (e.g., 
Polisensky and Ricotti, 2011), and as already mentioned WDM does not 


appear to be consistent with observed systematics of small galaxies. 

Another possibility is that the dark matter particles interact with them¬ 


selves much more strongly than they interact with ordinary matter (Spergel 


and Steinhardt, [2000 ). There are strong constraints on such self-interacting 


dark matter (SIDM) from colliding galaxy clusters (Harvey et ah, 2015 


Massey et ah, 2015), and in hydrodynamic simulations of dwarf galaxies 


SIDM has similar central cusps to CDM ( ]Bastidas Fry et 2015). SIDM 
can be velocity-dependent, at the cost of adding additional parameters, and 
if the self-interaction grows with an inverse power of velocity the effects can 


be strong in dwarf galaxies (Elbert et ah, 2014). An Aquarius-type sim¬ 


ulation but with velocity-dependent SIDM produced subhalos with inner 
density structure that may be compatible with the bright dSph satellites of 
the Milky Way (Vogelsberger et ah, 2012). Whether higher-resolution simu¬ 
lations of this type will turn out to be consistent with observations remains 
to be seen. 


1.6 Smaller Scale Issues: Dark Matter Halo Substructure 


The first strong indication of galaxy dark matter halo substructure was 
the flux ratio anomalies seen in quadruply imaged radio quasars (“radio 


quads”) ( 

Metcalf and Madau 

2001 

Dalai and Kochanek 

2002 

Metcalf and 

Zhao 2002 

). Smooth mass models of lensing galaxies can 

easily explain the 


observed positions of the images, but the predictions of such models of the 
corresponding fluxes are frequently observed to be strongly violated. Optical 
and X-ray quasars have such small angular sizes that the observed optical 
and X-ray flux anomalies can be caused by stars (“microlensing”), which 
has allowed a measurement of the stellar mass along the lines of sight in 
lensing galaxies (Pooley et ah, 2012). But because the quasar radio-emitting 


region is larger, the observed radio flux anomalies can only be caused by 
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relatively massive objects, with masses of order 10® to 10®Mq along the line 
of sight. After some controversy regarding whether ACDM simulations pre¬ 
dict enough dark matter substructure to account for the observations, the 
latest papers concur that the observations are consistent with standard the¬ 


ory, taking into account uncertainty in lens system ellipticity (Metcalf and 


Amara 2012) and intervening objects along the line of sight (Xu et al. 


2012 


2015). But this analysis is based on a relatively small number of observed 
systems (Table 2 of Chen et al. (2011) lists the 10 quads that have been 
observed in the radio or mid-IR), and further observational and theoretical 
work would be very helpful. 

Another gravitational lensing indication of dark matter halo substructure 
consistent with ACDM simulations comes from detailed analysis of galaxy- 


galaxy lensing (Vegetti et al., 2010, 2012, 2014), although much more such 


data will need to be analyzed to get strong constraints. Other gravitational 
lensing observations including time delays can probe the structure of dark 


Keeton and Moustakas 

2009 

)• 

Hezaveh et al. 


(2013, 2014) show that dark matter substructure can be detected using 


spatially resolved spectroscopy of gravitationally lensed dusty galaxies ob¬ 
served with ALMA. Nierenberg et al. (2014) demonstrates that subhalos can 
be detected using strongly lensed narrow-line quasar emission, as originally 
proposed by Moustakas and Metcalf (2003). 

The great thing about gravitational lensing is that it directly measures 
mass along the line of sight. This can provide important information that 
is difficult to obtain in other ways. For example, the absence of anomalous 
skewness in the distribution of high redshift Type la supernovae brightnesses 
compared with low redshift ones implies that massive compact halo objects 
(MACHOs) in the enormous mass range 10“^ to 10 ^®Mq cannot be the main 
constituent of dark matter in the universe (Metcalf and Silk, 2007). The low 
observed rate of gravitational microlensing of stars in the Large and Small 
Magellanic clouds by foreground compact objects implies that MACHOs in 
the mass range between 0.6 x 10“'^ and 15Mq cannot be a significant fraction 


of the dark matter in the halo of the Milky Way (Tisserand et ah, 2007). 


Gravitational microlensing could even detect free-floating planets down to 
10“®Mq, just one percent of the mass of the earth (Strigari et ah, 2012). 

A completely independent way of determining the amount of dark mat¬ 
ter halo substructure is to look carefully at the structure of dynamically 
cold stellar streams. Such streams come from the tidal disruption of small 
satellite galaxies or globular clusters. In numerical simulations, the streams 
suffer many tens of impacts from encounters with dark matter substruc¬ 
tures of mass 10® to IO^Mq during their lifetimes, which create fluctuations 
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in the stream surface density on scales of a few degrees or less. The observed 


streams contain just such fluctuations ( 

Yoon et al. 

Carlberg et al. 

2012 

Carlberg and Grillmair 

2013), 


2011; Carlberg, 2012 


evidence that the predicted population of subhalos is present in the halos 
of galaxies like the Milky Way and M31. Comparing additional observations 
of dynamically cold stellar streams with fully self-consistent simulations will 
give more detailed information about the substructure population. The Gaia 
spacecraft’s measurements of the positions and motions of vast numbers of 
Milky Way stars will be helpful in quantifying the nature of dark matter 
substructure (Ngan and Carlberg, 2014 Feldmann and Spolyar, 2015). 


1.7 Conclusions 

ACDM appears to be extremely successful in predicting the cosmic mi¬ 
crowave background and large-scale structure, including the observed distri¬ 
bution of galaxies both nearby and at high redshift. It has therefore become 
the standard cosmological framework within which to understand cosmo¬ 
logical structure formation, and it continues to teach us about galaxy for¬ 
mation and evolution. For example, I used to think that galaxies are pretty 
smooth, that they generally grow in size as they evolve, and that they are 
a combination of disks and spheroids. But as I discussed in Section 3, HST 
observations combined with high-resolution hydrodynamic simulations are 
showing that most star-forming galaxies are very clumpy; that galaxies often 
undergo compaction, which reduces their radius and greatly increases their 
central density; and that most lower-mass galaxies are not spheroids or disks 
but are instead elongated when their centers are dominated by dark matter. 

ACDM faces challenges on smaller scales. Although starbursts can rapidly 
drive gas out of the central regions of galaxies and thereby reduce the central 
dark matter density, it remains to be seen whether this and/or other baryonic 
physics can explain the observed rotation curves of the entire population 
of dwarf and low surface brightness (LSB) galaxies. If not, perhaps more 
complicated physics such as self-interacting dark matter may be needed. 
But standard ACDM appears to be successful in predicting the dark matter 
halo substructure that is now observed via gravitational lensing and stellar 
streams, and any alternative theory must do at least as well. 
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